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Abstract 

Popular density functionals for the exchange-correlation energy typically fail to reproduce the 
degeneracy of different ground states of open-shell atoms. As a remedy functionals which explic- 
itly depend on the current density have been suggested. We present an analysis of this problem 
by investigating functionals that explicitly depend on the Kohn-Sham orbitals. Going beyond the 
exact-exchange approximation by adding correlation in the form of the Colle-Salvetti functional 
we show how current-dependent terms enter the Colle-Salvetti expression and their relevance is 
evaluated. A very good description of the degeneracy of ground-states for atoms of the first and 
second row of the periodic table is obtained. 
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I. INTRODUCTION 



Common approximations to the exchange-correlation functional of density functional 
theory (DFT) jl[ 2] and spin-DFT (SDFT) y] often fail to reproduce the degeneracy of dif- 
ferent ground states. An illustrative example are ground states of open-shell atoms where 
one usually and erroneously obtains different total energies for states with zero and non- 
vanishing current density. The local spin density approximation (LSD A) gives rather 
small splittings (of the order of 1 kcal/ mole), but generalized gradient approximations 
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no. 



(GGAs) and meta-GGAs can introduce splittings of 10 kcal/mol |4,|s, 

These spurious energy splittings would vanish if the exact exchange-correlation func- 
tional could be used. The exchange-correlation energy functional can be represented 
in terms of the exchange-correlation hole function. Considering current-carrying states, 
Dobson showed how the expression for the exchange-hole curvature has to be changed 
by including current-dependent terms |12f | - Later, Becke found the same kind of terms 
in the short-range behavior of the exchange-correlation hole, and observed that they 
also enter the spin-like correlation-hole function |13Q . For open-shell atoms, inclusion 
of these current-dependent terms results in spurious energy splittings of less than 1 
kcal/mole [9]. Along these lines, Maximoff at al. [10] worked out a correction for the 
system-averaged exchange hole of the Perdew-Burke-Ernzerhof (PBE) GGA [14], which 
improves the corresponding spurious splittings. Alternatively, Tao and Perdew |ll[ 15 ] 
proposed a scheme for the extension of existing functionals using ideas of current density 
functional theory (CDFT) |16|, ll7l] which, again, improves the description of the degener- 
acy. 

The performance of the exact-exchange (EXX) energy functional - which, by definition, 
describes the exchange hole correctly - has been evaluated for the spurious splittings in 
DFT and SDFT [18]. In the EXX-DFT (i.e., spin-restricted calculations using one and the 
same Kohn-Sham potential for spin-up and spin-down orbitals) the degeneracy is well re- 
produced to within 0.6 kcal/mole but, surprisingly, in EXX-SDFT (i.e., spin-unrestricted 
calculations using two Kohn-Sham potentials, one for spin-up and one for spin-down or- 
bitals) spurious splittings up to 3 kcal/mole are obtained. In particular, current-carrying 
states always have higher total energies than states without current. This observation 

nnn, 

motivated the applications of the optimized-effective-potential (OEP) method [19, 201,1211] 
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generalized to current-spin-density functional theory (CSDFT) to these current-carrying 
states I22I] . As expected, EXX-CSDFT total energies for current-carrying states are lower 
than those of EXX-SDFT. However, this lowering is too small to give a substantial im- 
provement of the spurious energy splittings. These studies lead to the conclusion that 
correlation is needed for any further improvement. 

The construction of a correlation energy functional compatible with EXX is a difficult 
task 122, 



24, 



250, but for spherical atoms it was found that EXX combined with the Colle- 



Salvetti (CS) functional for correlation 



27, 



28(] leads to very accurate total energies 



291 ] . The CS functional has been used to derive the popular Lee-Yang-Parr (LYP) func- 
tional 1 3(3], which is most commonly used together with Becke's exchange functional [31] 



32, 



(BLYP) and in hybrid schemes such as B3LYP 
tion energy functional also has its limitations [34, 
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. On the other hand, the CS correla- 



36Q . In particular, while short-range 



correlations are well described |35Q , very important long-range correlations are missing. 
These correlations often cannot be ignored in molecules and solids, but are negligible in 
atoms. This fact, together with the encouraging results for spherical atoms [Q], indicates 
that it is appropriate to employ the CS functional to analyze the degeneracy problem for 
open-shell atoms beyond EXX. Furthermore, the expression of the CS functional also al- 
lows a reconsideration of the relevance of the orbital currents as ingredient of correlation 
functionals. 

Although the general density functional formalism to deal with degenerate ground 
states includes densities which can only be obtained by a weighted sum of several deter- 
minantal densities 1 37, 38], as in many previous investigations [gll^[l^[lij] we only con- 
sider densities which may be represented by a single Slater determinant of Kohn-Sham 
orbitals. 



II. THEORY 



Going beyond the EXX approximation, we here consider the correlation-energy func- 



tional of Colle and Salvetti |26Q . This expression relies on the assumption that the corre- 
lated two-body reduced density matrix may be approximated by the Hartree-Fock (HF) 
two-body reduced density matrix pf F (ri,r 2 ), multiplied by a Jastrow-type correlation 
factor. After a series of approximations, the following expression is obtained for the cor- 
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relation energy 



E r = -4a / dr 



pf F (r,r) 


l + tyr§(r) 


(r+f ,r-l ) 


s=0- 


e -cp 3(r) 


P(r) " 


l + dp-s(r) 



(1) 



where p^ F (r, s) is expressed in terms of the average and relative coordinates r = |(rx+r 2 ) 
and s = ri — r 2 . Here, p(r) is the electron density and the constants a = 0.049, b = 0.132, 
c = 0.2533, d = 0.349 are determined by a fitting procedure using the Hartree-Fock (HF) 
orbitals for the Helium atom. 

Following Lee, Yang and Parr, this expression can be restated as a formula involving 
only the total charge-density, the charge-density of each Hartree-Fock orbital and their 
gradient and Laplacian ||30j]. In this derivation, the single-particle orbitals are tacitly as- 
sumed to be real. We denote the resulting expression as CSLYP. In the following, we 
relax this restriction and consider complex orbitals. We then proceed in analogy to the 
inclusion of current-dependent terms in the Fermi-hole curvature [9, 



121 ] and in the ex- 



tension of the electron-localization-function (ELF) |39fl for time-dependent states [40]. As 
a consequence, in addition to the term already present in CSLYP expression, the current 
densities of the single-particle orbitals appear in the final formula. In order to obtain this 
expression, which in the following will be denoted as JCSLYP, we rewrite the Laplacian of 
the Hartree-Fock (HF) two-body reduced density matrix in Eq.([l]) in terms of the original 
particle coordinates 



s s 
r H — , r 

2 ' 2 



~V? + jV 2 2 -^V 1 -V 2 )pf ir (r J .r 2 )| ) ,. 1 ., 



where 



Here, 



pHn,r 2 ) = kr0p(r 2 ) -^^(l.'JCK^) ■ 



P?f(ri,r 2 ) = ^^(ri)^, ff (r 2 ) 



(2) 



(3) 



(4) 



fc=i 



is the first-order HF density matrix (for a single Slater determinant) expressed in terms 
of the single-particle orbitals i^k,a{ r )- The corresponding spin-density is simply given by 



Ar(r) = p ha (r,r 



(5) 
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Allowing the single-particle orbitals iftk,a( T ) to be complex, a given orbital not only gives 
the contribution Pfc )(T (r) = |^ fc (7 (r)| 2 to the density, but also the contribution j p k,a( r ) — 
Im (ip k>a (r)Vipl CT (r)) to the paramagnetic current density which is given by 



(6) 



k=l 



After some straightforward algebra, the Laplacian of the second-order HF reduced den- 
sity matrix takes the final form 



2 JUF 



r H — , r 

2 ' 2 



i p(r) V 2 p(r) - ~ (Vp(r)) 2 - ~ £ p CT (r) V 2 p CT (r) (7) 
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^ (Vp fej(T (r)) 2 



fc=i 



where 



J(r) = ^p CT (r) 



i 2 ( 

J per V 



Na i 2 (VI 



(8) 



contains all the current-dependent terms. Alternatively, Eq. (0 may also be expressed in 
terms of the non-interacting kinetic energy density 



as 



2 JUF 



k=l 



r H — , r 

2 ' 2 



fe=l 
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Pk,a( r ) 



, AT CT 



2 fcj P*» 



(9) 



s=0 



i p(r) V 2 p(r) - ~ (Vp(r)) 2 - ~ £ p CT (r) V 2 p CT (r) 

cr 

+ ^(2 P(J (r)r CT (r)-j p(7 (r) 2 ). (10) 



Comparison of Eq.© and Eq. ((T0)) shows that J(r), as defined in Eq.©, also contains 
current-dependent terms coming from the kinetic energy density. Thus, this would also 
suggest to reconsider the gradient expansion of r for current-carrying states. Impor- 
tant consequences maybe expected for all approximated exchange-correlation functional 
involving the kinetic energy density as ingredient, such as the CS functional, and meta- 
GGAs. This issue will be specifically considered in a future work. In the next section, we 
assess the performance of the CS functional, and, in particular, the relevance of J(r), in 
reproducing the degeneracy of atomic states. 
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III. RESULTS AND DISCUSSION 



We consider ground states of open-shell atoms having densities that can be repre- 
sented by a single Slater determinant of Kohn-Sham orbitals. Due to the symmetry of the 
problem, these Slater determinants are eigenstates of the z-component of both spin and 
orbital angular momentum. The single-particle orbitals are generally complex-valued 
and states with different total magnetic quantum numbers, M L , correspond to differ- 
ent current densities. By means of an accurate exchange-correlation functional the same 
total energies would be obtained. In the previous section, we have shown how current- 
dependent terms enter the expression of the CS functional when complex-valued orbitals 
are considered. Here, we evaluate the performance of EXX plus the CS functional in re- 
producing the degeneracy and study the effect of the current-dependent terms, in SDFT 
and DFT calculations. 

We consider atoms of the first and second row of the periodic table: these are the ref- 
erence cases for which a vast amount of numerical data is available 
In analogy to the procedure where Hartree-Fock orbitals are used as input to the CS for- 
mula, we have evaluated the correlation energies in a post-hoc fashion using Kohn-Sham 
(KS) orbitals. We expand the KS orbitals in Slater-type basis functions (QZ4P of Ref. ^42]) 
for the radial part, multiplied with spherical harmonics for the angular part. We obtain 
the KS orbitals from self -consistent EXX-only calculations employing the approximation 
of Krieger, Li, and Iafrate (KLI) [41|], which has been shown to be extremely good at least 
for small systems [21]. In principle a functional should be evaluated with KS orbitals 
obtained from self-consistent calculations, and this is certainly possible for the CS func- 



tional |29|] . However, thanks to the variational nature of DFT, it is common experience to 



observe only minor quantitative differences between post-hoc and self-consistent evalu- 



for solving the 



ation of total energies. This is the reason why the functionals designee 
degeneracy problem are typically evaluated in a post-hoc manner |J9|, llOi llll ] 

Table H] shows the spurious energy splittings (difference in the total energies) between 
Kohn-Sham Slater determinants with total magnetic quantum number \M L \ = 1 and 
M L = 0, from our SDFT and DFT calculations. The deviation of these total energies from 
the exact values is plotted in Fig. (1), where again the spurious energy splittings are visi- 
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SDFT (DFT) 


Atom 


&-JCSLYP 


&-CSLYP 


B 


0.8 (-0.3) 


2.4 (1.4) 


C 


0.9 (-0.1) 


-3.2 (-4.3) 


O 


-0.6 (-1.9) 


0.9 (-0.4) 


F 


-0.1 (-1.5) 


-3.5 (-5.1) 


Al 


0.4 (-0.5) 


1.1 (0.2) 


Si 


0.5 (-0.4) 


-1.2 (-2.2) 


S 


0.1 (-1.6) 


1.1 (-0.7) 


CI 


0.7 (-1.3) 


-1.1 (-3.2) 


me 


0.3 (-1.0) 


-0.4 (-1.8) 


mae 


0.5 (1.0) 


1.8 (2.2) 



TABLE I: Spurios energy splittings, A = E(\Ml\ = 1) — E(Ml = 0) in kcal/mol for open-shell 
atoms, computed in SDFT (DFT results in parenthesis for comparison). Correlation energy has 
been added to the KLI-EXX energies including (JCSLYP) and neglecting (CSLYP) the current terms 
of Eq. 10. The last row shows the mean error (me) and mean absolute (mae) of the spurious 
splittings. 



ble. These results highlight the importance of including J(r) in Eq. 0. In particular, it is 
remarkable to observe that SDFT splittings are within 0.9 kcal/ mol (with a mean error of 
0.3 kcal/ mol and a mean absolute error of 0.5 kcal/ mol), and the corresponding DFT spu- 
rious energy splittings are less than 1.9 kcal/ mol (with mean errors of 1.0 kcal/ mol). It is 
worthwhile to note that in several cases inclusion of correlation leads to current-carrying 
states (| Ml | = 1) with lower total energy than zero-current states (M L = 0). These results 
are in contrast to EXX-only Q18Q cases where: (a) the zero-current states are always lowest 
in energy and (b) the spurious energy splittings are always smaller in DFT than in SDFT. 

Going beyond EXX by including correlation in the form of CS functional accurate to- 
tal energies can be obtained within the OEP method uM- Figure dD and Table [U show 
the deviations from exact total energies for the states with different magnetic quantum 
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FIG. 1: Deviation from exact total energies for SDFT and DFT calculations employing the CS 
functional, including (JCSLYP) and not including (CSLYP) the current-dependent term J in Eq. (0. 
States with different magnetic quantum numbers Mr are plotted. Exact total energies are taken 

n 

from Ref. [29] and references therein. 

numbers, i.e. different current-carrying states. This further emphasizes the importance 
of proper inclusion of J(r) in Eq. Q. 



IV. CONCLUSIONS 



We have shown that going beyond the exact-exchange approximation by including 
correlation energy in the form of the Colle-Salvetti functional leads to a very good de- 
scription of the degeneracy of open-shell atoms in both SDFT and DFT calculations. 
Comparing DFT and SDFT results for the first and second row of the periodic table, on 
average we observe a reduction of the spurious energy splittings and better total ener- 
gies for SDFT. Furthermore, we have also shown how current-dependent terms enter the 
expression of the Colle-Salvetti functional. If these terms are neglected, the degeneracy 
is not well described and the total energies are also less accurate. Thus, this analysis re- 
confirms the advantage of properly including the orbital current dependent terms as an 
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JCSLYP (CSLYP) 




SDFT 


DFT 


SDFT 


DFT 


\Ml\ 








1 


1 


me 


2.3 (16.1) 


4.7 (18.5) 


2.6 (15.7) 


3.7 (16.8) 


mae 


4.8 (16.4) 


5.4 (18.5) 


4.7 (15.7) 


4.7 (16.8) 



TABLE II: Mean error (me) and mean absolute error (mae) in the total energies for the CS 
functional, including (JCSLYP) and not including (CSLYP results in parenthesis) the current- 
dependent term J in Eq. (fT)), in kcal/ mol. Exact total energies are taken from Ref. \z9] and ref- 
erences therein. 

ingredient in correlation functionals. 
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